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Proofs that Fredholm determinants of transfer operators for hyperbolic flows are entire can be ex- 
tended to a large new class of multiplicative evolution operators. We construct such operators both 
for the Gutzwiller semi-classical quantum mechanics and for classical thermodynamic formalism, 
and introduce a new functional determinant which is expected to be entire for Axiom A flows, and 
whose zeros coincide with the zeros of the Gutzwiller- Voros zeta function. 



It has been established recently [yy,|j that the 
Gutzwiller- Voros zeta function derived from the 
Gutzwiller semi-classical trace formula is not an en- 
tire function. In this letter we construct a new classical 
evolution operator whose Fredholm determinant ( p^ ) 
is entire for Axiom A flows, and whose spectrum contains 
the semi-classical Gutzwiller spectrum. Many physically 
realistic chaotic scattering systems [^,0 are of Axiom A 
type, and for them the new determinant has much bet- 
ter convergence properties than the Gutzwiller- Voros and 
Ruelle type zeta functions utilized previously [@,||. 

The main idea, extending the dynamical system to the 
tangent space of the flow, is suggested by one of the 
standard numerical methods for evaluation of Lyapunov 
exponents start at xq with an initial vector ^(0), 
and let the flow transport it to ^(t) along the trajectory 
x{t) = /*(xo). The growth rate of this vector is multi- 
plicative along the trajectory 
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and can be represented by the trajectory of a "unit" vec- 
tor u(t) multiplied by the factor For asymp- 
totic times and for almost every initial (xq, this fac- 
tor converges to the leading eigenvalue of the linearized 
stability matrix for the flow. 

We implement this multiplicative evaluation of stabil- 
ity eigenvalues by adjoining [l0| the d-dimensional trans- 
verse tangent space ^ G TUx^ £,{x) ■ 'v{x) = 0, to the 
((i-l-l)-dimensional dynamical evolution space x € U C 
rpj-^g dynamics in the (a;,i^) G [/ x TUx space is 
governed by the system of equations of variations [0 : 

X = v(x) , ^ = Dv(a;)^ . 

Here Dv(a;) is the transverse derivative matrix of the 
flow. We write the solution as 



x{t) = f{xo), m^^'M-^o, 



(2) 



with the tangent space vector ^ transported by the trans- 
verse stability matrix J* (xq) — dx{t) / dxo . In order to de- 
termine the length of the vector ^ we introduce a signed 
norm, a differentiable scalar function g{^) with the prop- 
erty .g(A^) — A(7(^) for any number A. While in general 
such norm is a space dependent function g{^, x), we shall 
assume here for reasons of notational simplicity that g is 
a function of ^ only. An example is the function 
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Any vector ^ e TUx can now be represented by the prod- 
uct ^ = Au, where u is a "unit" vector in the sense that 
its signed norm is ^(u) = 1, and the factor 

A\xo, uo) = gim) = .9(J*(xo) • uq) (4) 

is the multiplicative "stretching" factor introduced in (|l|) 

A*'+*(xo,uo) = A''{x{t),u{t))A\xo,uo). 

The u evolution constrained to ETg^x, the space of unit 
tangent vectors transverse to the flow v, is given by 
rescaling of (0): 



u' = R\x,\i) 
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A*(a;, u) 



J*(a;) -u. 



(5) 



Eqs. (0), and (g|) enable us to define a multiplicative 
evolution operator on the extended space U x ETg x 



C\u', x'; u, x) = e''("'<5(2:' - f{x)) 



^5[m' ~ R\x,u)) 
' |A*(a;,u)|/3-i 



(6) 



where /i is a function additive along the trajectory, and 
/3 is a number. This should be contrasted to "thermody- 
namic" operators of form 



1 



|A*(x)|/3- 



with A*(a;) an eigenvalue of J*(x). Such operators are 
not multiplicative in two or more transverse dimensions, 
for the simple reason that the eigenvalues of successive 
stability matrices are in general not multiplicative 

Aab ^ AaAb- 

Here Job = JaJb is the stability matrix of the trajectory 
consisting of consecutive segments a and 6, 3a and J;, are 
the stability matrices for these segments separately, and 
A's are their eigenvalues. In particular, this lack of mul- 
tiplicative property for A's had until now frustrated at- 
tempts to construct evolution operators whose spec- 
trum contains the semi-classical Gutzwiller spectrum. 

In order to derive the trace formula for the operator (||) 
we need to evaluate Tr£* = J dxduC'^{u,x;u,x). The 
J dx integral yields [ p^ a weighted sum over primitive 
periodic orbits p and their repetitions r: 



Tr £* = ^ Tp 
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du 



^ I det(l - J;) I 

5{\i- R^p'^ixp.n)) 
|A^p'-(.Tp,u)|/3-i ' 
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where Jp is the prime cycle p transverse stability matrix. 
As we shall see below, Ap^^ is intrinsic to cycle p, and 
independent of any particular cycle point Xp. 

We note next that if the trajectory is periodic 

with period T, the tangent space contains d periodic so- 
lutions 



e^(x(T + t)) = e,(a;(t)), i = 1, 



corresponding to the d unit eigenvectors {ei, 62, • • • , e^;} 
of the transverse stability matrix, with "stretching" fac- 
tors given by its eigenvalues 

Jp(x) • e^{x) = Ap^j ei{x) , z = 1, ...,d. 

The / du integral in picks up contributions from these 
periodic solutions. In order to compute the stability of 
the j-th eigendirection solution, it is convenient to expand 
the variation around the eigenvector in the stability 
matrix eigenbasis = ^ 5u^ ei . The variation of the 
map (0) at a complete period t — T is then given by 



J ■ (5u J ■ 

g{J ■ e.i) g{3 ■ 6^)2 
A, 
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p.k f _ ^ dgje. 



du 

Suk ■ 



J Su 
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The 5ui component does not contribute to this sum since 
g{ei + duiGi) = 1 -|- dui implies dg{ei)/dui — 1. Indeed, 
infinitesimal variations du must satisfy 



g{u + du) = g{u) = 1 



due 



so the allowed variations are of form 



<5U = ^ ( Gfc 

k^i 



dg{e,) 
ouk 



\Ck\ < 1, 



and in the neighborhood of the eigenvector the J u 
integral can be expressed as 



du = 



Yl^dck 

k^i 



Inserting these variations into the / du integral we obtain 

/ dud{e, +Su- R^{e,) - 6R^{e,) + . . .) 
•'g 

- [Y[dckS ((1 - Afc/A,)cfc + . . .) 

k^i 

^n^^AiTAd' 

and the / du trace ^ becomes 

A„ _ = V I TT I 

^|A:,|/5-ili|i_A;,/A^J 



(9) 



k=^i 



The corresponding Fredholm determinant is obtained by 
observing |16| that the Laplace transform of the trace 



TtC{s) = / 
J04 



dt e^' Tr C{t) 



is a logarithmic derivative Tr£(s) = — ^ logF(s) of the 
Fredholm determinant: 



f^(/3,s)=expU^ 



4hp+sTp)r 



r I det(l - J;) 



rAp,.(/3) . (10) 



This determinant is the central result of this paper. Its 
zeros correspond to the eigenvalues of the evolution op- 
erator (^), and can be evaluated by standard cycle ex- 
pansion methods [p[p^]. 

In the "thermodynamic" formalism for clas- 

sical chaotic systems, and in the Gutzwiller semi-classical 
description of systems with chaotic classical counterpart 

, /3 is a parameter which plays the role of "inverse tem- 
perature" , and hp is the integral of some weight function 
h{x) evaluated along the prime periodic orbit p. In the 
semi-classical quantization case |5| 



/3=l/2, 



iSp/h + WpTT/2, 



0, 



(11) 



where Sp is the action of the periodic orbit, and Vp its 
Maslov index. The classical correlation spectra are given 
by /3 = 1 and hp = 0. 
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The simplest application of (|T^) is to 3-dimensional hy- 
perbolic Hamiltonian flows (higher dimensions need fur- 
ther evolution operators, for outer products of vectors 
rather than single vectors). In this case Ai = I/A2 = A, 
and the Fredholm determinant is given by 



F„{P,s) = exp 



cr. 



P 



„r(/ip+sTp) 



Ap,.(/3) 



iA-i-^+^ iA-r3 

1 - l/A^r + 1 - 1/A2r 



(12) 



The extra multiplicative factor is set to the eigenvalue 
sign a — A/|A| for F_, and to = 1 for this will be 
used below. 

The Gutzwiller-Voros zeta function corresponds to set- 
ting Ap,r = |Ap|^/2(l - 1/Ap, and the "quantum Fred- 
holm determinant" |^ is obtained by setting Ap^r = 
|A^|i/2 in (^. The practical advantage of (0) over 
the more familiar Gutzwiller-Voros and Ruelle type zeta 
functions was demonstrated by detailed numerical stud- 
ies Q of the related quantum Fredholm determinant . 
In the systems studied, the quantum Fredholm determi- 
nant appeared to be entire for all practical purposes; only 
the most recent numerical investigations [T^ reveal poles 
absent in (p^, see fig. l} 

It can be shown |l| ] that the Fredholm determinant 
obtained by keeping only one of the terms in the sum in 
(^ is entire. This enables us to show that the Gutzwiller- 
Voros zeta function Z{E) for Axiom A flows is meromor- 
phic in the complex E plane, as it can be written as a 
ratio of entire functions; for 2-dimcnsional Hamiltonian 
systems 



F+{\,E)F^{l,E) 
Fa,E)F+[lE)' 



(13) 



where F„ includes only the first term in the Ap^^ sum 
(p^). The zeros of the Gutzwiller-Voros zeta function 
coincide with the ones obtained from F+(i,i?). Such 
relations follow by inserting into Ap^^ identities like 



1-1/A'- A'-(l-l/A'^) 

(ct weight in ( |T^ ) is needed to account for the 1/A = <t/\A\ 
term in the above indentity). 

We illustrate a choice of the g{^) function and construc- 
tion of the R dynamics (^) by a simple explicit exam- 
ple: 2-dimensional Hamiltonian dynamics reduced to a 
2-dimensional Poincare section return map Xi+i = f{xi). 
The stability matrix of cycle p is a product of the [2 x 2] 
stability matrices 

where Aj = dfi{xi)/dxi, and so on. Assume the signed 
norm (0) and multiply an initial unit vector by the first 



stability matrix in the product. The resulting vector can 
be written as 



1 



= (CiKi+Di) ( 



Hence the dynamics acts on the unit vectors as a rational 
fraction transformation 



Kk+l = R{xk, Kk) 



AkKk + Bk 



(14) 



with the signed norm (^) of the iterated vector given by 

Up 

A^-{xi,Kl)^\[{C,K, + Di). 

i=l 

(In the case of 2-dimensional billiards, k is the 
Bunimovich-Sinai curvature [p9[). For a periodic or- 
bit Knp — kq, the unit vector is an eigenvector of the 
stability matrix, and the corresponding eigenvalue is 
Ap = n(Ci'^i + ^i)- For 3-dimensional hyperbohc flows 
there are two distinct Kp^i values, one for the expand- 
ing, and the other for the contracting eigenvalue. The 
derivative of i?"" is easily evaluated; 
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—R^^{Ki,xi) = l[ R'{n,,Xi) = W 



1 



(we have used dct(Ji) = 1) and we again obtain the trace 
formula (|l|). 

In conclusion, we have constructed a classical evo- 
lution operator for semi-classical quantization, and de- 
rived a new determinant for periodic orbit quantization 
of chaotic dynamical systems. The main virtue of the de- 
terminant ( p^ ) is that the theorem of H.H. Rugh ||2^, ap- 
plicable to multiplicative evolution operators such as (|^), 
implies that this determinant should be entire for Axiom 
A flows, ie. free of poles in the entire complex s or com- 
plex energy plane. One consequence of this general result 
is that the cycle expansions of the new Fredholm deter- 
minant should converge faster than exponentially. Our 
numerical tests support above claims [|l8|; one example 
is given in fig. || which demonstrates that the new de- 
terminant is analytic to the limit of numerical precision 
of current cycle expansions, and well beyond both the 
Gutzwiller-Voros zeta function and the quantum Fred- 
holm determinant. 
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FIG. 1. The leading semi-classical resonances in the k com- 
plex wave- number plane (x) for the determinant (|l^, com- 
pared with (□) the zeros of the "quantum Fredholm deter- 
minant" I^Q. While the quantum Fredholm determinant 
has a finite region of analyticity (the bottom line of zeros 
reflects a pole expected at Im{k) = —4.559843..., indi- 
cated by the dashed line labeled "QF"), the new determi- 
nant shows no numerical indication of any poles, and en- 
ables us to reach resonances deeper down in the complex 
plane. The Gutzwiller- Voros determinant is reliable only 
down to the line labeled "GV", Im(fc) = -2.491905..., the 
upper bound on the poles of determinant _F_(|,_E) in (p^. 
The dynamical system tested is the Hamiltonian Henon map 
Xk+i = 1 — ax^. — Xk-i at a — 20, with cycle expansions 
truncated to cycles up to period 18. We take as action 
Sp/h — kup, and as the Maslov phase z/p — 2np (this is a 
normal-form model of a 3-disc repeller, see ref. 
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